Casimir amplitudes and capillary condensation of near-critical fluids between parallel 

plates: Renormalized local functional theory 
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We investigate the critical behavior of a near-critical fluid confined between two parallel plates in 
contact with a reservoir by calculating the order parameter profile and the Casimir amplitudes (for 
the force density and for the grand potential). Our results are applicable to one-component fluids 
and binary mixtures. We assume that the walls absorb one of the fluid components selectively for 
binary mixtures. We propose a renormalized local functional theory accounting for the fluctuation 
effects. Analysis is performed in the plane of the temperature T and the order parameter in the 
reservoir ipoa- Our theory is universal if the physical quantities are scaled appropriately. If the 
component favored by the walls is slightly poor in the reservoir, there appears a line of first-order 
phase transition of capillary condensation outside the bulk coexistence curve. The excess adsorp- 
tion changes discontinuously between condensed and noncondensed states at the transition. With 
increasing T, the transition line ends at a capillary critical point T = T^^ slightly lower than the 
bulk critical temperature Tc- The Casimir amplitudes are larger than their critical-point values by 
10-100 times at off-critical compositions near the capillary condensation line. 



I. INTRODUCTION 

Understanding the phase behavior of fluids confined in 
narrow regions is crucial in the physics of fluids in porous 
media and for surface force experiments^. It strongly 
depends on the molecular interactions. In binary mix- 
tures, one of the solvent components is preferentially at- 
tracted to the wall^i^. In liquid water in contact with 
a hydrophobic surface, water molecules tend to be sep- 
arated from the surface due to the hydrogen bonding 
among the water molecules, resulting in the formation 
of a gaseous layer on the surface^. In aqueous mixtures 
with salt, prewetting behavior is much intensified by the 
composition-dependent surface ionization and the pref- 
erential solvation of ions^"— . In these examples, hetero- 
geneities in the density or the composition are induced 
near the wall, often resulting in wetting or drying transi- 
tion on the wall. The confinement effects can be dramatic 
when the spatial scale of the fluid heterogeneities is on 
the order of the wall separation. 

Narrow regions may be filled with the phase favored 
by the confining walls or may hold some fraction of the 
disfavored phaseii^iiSrii^, depending on the temperature 
T, and the reservoir chemical potential difference ^ocr^ 
(or the reservoir order parameter ■f/'oo) for each given wall 
separation D between these two states, there is a first- 
order phase transition line slightly outside the bulk co- 
existence curve in the T-fj,oc plane. In this paper, we 
call it a capillary condensation line, though this name is 
usually used for gas- liquid phase separation in narrow re- 
gions. For binary mixtures, there is a discontinuity in the 
preferential adsorption and the osmotic pressure across 
this line. Maciolek et al.^^ numerically found a capillary 
condensation line in two-dimensional Ising films. Samin 
and Tsori^^ found a first-order transition in binary mix- 
tures with ions between parallel plates including the se- 
lective solvation in the mean-field theory, where the os- 
motic pressure is discontinuous. Our aim in this paper 



is to investigate this capillary condensation transition in 
near-critical fluids between parallel plates, taking into 
account the renormalization effect of the critical fluctua- 
tions. For simplicity, no ions will be assumed. 

In near-critical fluids, the wall-induced heterogeneities 
extend over mesoscopic length scales, leading to an in- 
triguing interplay between the finite-size effect and the 
molecular interactionsi^ii^. Note that various surface 
phase transitions are known near the bulk criticality de- 
pending on the type of surfaceii"— . According to the 
prediction by Fisher and de Gennes, the free energy of a 
film with thickness D and area S at the bulk criticality 
consists of a bulk part proportional to the volume SD, 
wall contributions proportional to the area 5', and the 
interaction part 
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AF = -SknTrD- 



-d+l 



A, 



(1.1) 



where d is the space dimensionality, and A is a univer- 
sal number only depending on the type of the boundary 
conditions. In the literature, A has been called the the 
Casimir amplitude^^. Positivity (negativity) of A implies 
attractive (repulsive) interaction between the plates. The 
force density between the plates is then written as 



— AF = SksT^D-U, 
oD 



(1.2) 



where A is another amplitude. The amplitudes A and A 
depend on the reduced temperature r — T /Tc — 1 and 
the reservoir chemical potential difference /.too. However, 
A and A have been measured as functions of r along 
the critical path fXoo = O^i"— . In this paper, we cal- 
culate them in the r-fiao plane to find their dramatic 
enhancement near the capillary condensation line, where 
the reservoir composition is poor in the component fa- 
vored by the walls. In accord with our result, A increased 
dramatically around a capillary critical line under a mag- 
netic field in two-dimensional Ising filmsi^. 
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In binary mixtures, the adsorption-induced composi- 
tion disturbances produce an attractive interaction be- 
tween sohd objects^i^. Indeed, reversible aggregation of 
colloid particles have been observed close to the criti- 
cal point at off-critical compositions^"—. On approach- 
ing the solvent criticality, this interaction becomes long- 
ranged and universal in the strong adsorption limi t^^'^^ . 
Recently, the near-critical colloid-wall interaction has 
been measured directly at the critical composition^SiiL. 
However, the interaction should be much larger at off- 
critical compositions. We mention a microscopic theory 
by Hopkins et al^, who found enhancement of the col- 
loid interaction in one-phase environments poor in the 
component favorted by the colloid surface. Furthermore, 
ions can strongly affect the colloid interaction in an aque- 
ous mixture^, since the ion distributions become highly 
heterogeneous near the charged surface inducing forma- 
tion of a wetting layer on the surface. 

The organization of this paper is as follows. In Sec. II, 
we will present a local functional theory of a binary mix- 
ture at the critical temperature t = T/Tc — 1 = 0. The 
Casimir amplitudes can easily be calculated from this 
model. In Sec. HI, we will extend our theory in Sec. II to 
the case of nonvanishing reduced temperature r account- 
ing for the renormalization effect. We can then calculate 
the Casimir amplitudes for nonvanishing r and and 
predict the capillary condensation transition. 



II. CRITICAL BEHAVIOR AT T = 

In this section, we treat a near-critical binary mixture 
at the critical temperature T = Tc at a fixed pressure 
without ions. Using the free energy proposed by Fisher 
and Au-Yangi^ at T = Tc, we calculate the composition 
profiles and the Casimir amplitudes in the semi-infinite 
and film geometries. We mention similar calculations by 
Borjan and Upton^li^ii^ along the critical path /Xqo = 0. 
The method in this section will be generalized to the case 
of nonvanishing /ioo and T — Tc in the next section. 



A. Model of Fisher and Au-Yang 

The order parameter is the deviation tp = c — Cc oi the 
composition c from its critical value Cc- Supposing spatial 
variations of if) on scales much longer than the molecular 
diameter, we use the following form for the singular free 
energy density including the gradient contributioni^. 



/loc = /W + ^C(V)|V^|2, 

The total singular free energy is given by 



F 



(2.1) 



(2.2) 



where /Zqo is a given chemical potential difference in the 
reservoir— (magnetic field h for Ising spin systems) . The 



second term is the surface integral of the surface free 
energy fsi'4') arising from the short-range interaction be- 
tween the fluid and the wall^ (with dS being the surface 
element). At T = Tc, f{ip) and C{ip) are proportional to 
fractional powers of as 



cm 



(2.3) 
(2.4) 



where Bq and Cq are positive constants. The S, -q, v, and 
/3 are the usual critical exponents for Ising-like systems. 
For the spatial dimensionality d in the range 2 < c? < 4, 
they satisfy the exponent relations^-. 



d + 2~r] 
d- 2- 



dv 
J 



1, 



d-2 



V 



(2.5) 



The exponent rj is very small in the range 0.03 — 0.04 
for d — 3. The other exponenta we will use are 7 = 
(2 — 77)1/ = (6 — l)/3 and a = 2 — vd. In our numerical 
analysis to follow, we set a = 2 — vd ^ 0.11, /? — 0.325, 
5 = 4.815, V = 0.630, and 77 = 0.0317. They are three- 
dimensional values satisfying Eq.(2.5). 

We define the local correlation length ^ioc('0) by 



6oc(V') = {2C^lB^Y/^\^\-''/P. 



(2.6) 



The two terms in Eq.(2.1) are of the same order if we set 
IVtAP ~ V'V^ioc- The combination {2Cq/BqYI'^ is a mi- 
croscopic length if t/' denotes the composition deviation. 
In terms of ^loc, we rewrite Eq.(2.1) asi^ 

Cfoc/ioc/fcBTc = Ac[l + eL|V^|V4V'']. (2.7) 

Here the coefficient Ac is expressed as 

Ac = Bo{2Co/Bnf/^/kBTc. 

= {AC^/BokBTc)iBo/2CoY^\ (2.8) 

Hereafter e = 4 — d. In this model free energy, the 
critical fiuctuations with wave lengths shorter than ^loc 
have already been renormalized or coarse-grained. Thus 
minimizing F in Eq.(2.2) yields the average profile of 
near the walls, where we neglect the thermal fiuctuations 
with wavelengths shorter than ^loc. The theory of the 
two-scale- factor universality indicates that Ac in Eq.(2.8) 
should be a universal number (independent of the mix- 
ture species The second line of Eq.(2.8) indicates 
that its e expansion form is Ac — IS/n'^e + ■ ■ ■ as e 0. 
In the next section below Eq.(3.16), we will estimate Ac 
to be 1.49 at d = 3. 

For a thin film with thickness D, however, the local 
forms in Eqs.(2.3) and (2.4) are not justified in the spatial 
region where ^ioc(V') ^ D holds, since the length of the 
critical fiuctuations in the perpendicular direction cannot 
exceed D. Here we should note that two-dimensional 
critical fiuctuations varying in the lateral plane emerge 
near the critical point of capillary condensation, though 
they are neglected in this paper (see Sec. HI). From the 
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balance ^ioc(V'-d) = D, we may introduce a characteristic 
order parameter by 



(2.9) 



Then IV'I/V'd = (D / ^loc)'^ ^ " ■ In the film geometry ip will 
be measured in units of ipD- 



B. Profiles in tlie semi-infinite geometry 

We first consider the semi-infinite system (0 < z < cxi) 
to seek the one-dimensional profile ip = ^{z). Here, 
■0(0) = V'o > at z = and i^{z) ipao as 2 — 00. 
We treat ipQ and V'cx) as given parameters not explicitly 
imposing the usual boundary condition Cdip/dz = fi. 
at z = 0^, where = dfs/dip with fs being the sur- 
face free energy density in Eq.(2.2) (see comments below 
Eq.(2.17)). We shall see how the strong adsorption limit 
"00 ~^ 00 is attained near the criticality. 

The chemical potential far from the wall is written as 



(2.10) 



where the second line follows from Eq.(2.3). In equilib- 
rium, we minimize F to obtain 

/i(0) - c{i;)r - \c'm^'? = A^oo, (2.11) 

where V' = d^j/dz, iJj" = d^i^/dz^, and C'{ip) = dC/dip. 
This equation is integrated to give^ 



(2.12) 



Here LUsitp) is the excess grand potential density in the 
semi-infinite case defined by 

^ /(-ip) - foo - fJ-oo{'>Jj ~ Ipoo)- (2.13) 

where /cx) = fiipao)- We further integrate Eq.(2.12) as 

(2.14) 



lip 



where the upper bound is ipo = V'(O) and the lower bound 
tends to ipoc as z increases. 

In particular, for ipoo = or at the criticality in the 
bulk, we simply have uJs = f = Bo\'ip\^'^^ . We may read- 
ily solve Eq.(2.12) as^ 

i,iz)^i,o/[i + z/eof/'' 

- iP/2i.f/^2Co/Bof/^''iz + 4)-^/^ (2.15) 
where £q is the transition length related to tpo by 

£0 = iP/i^){Co/2Boy^'iPr^^- (2-16) 

From Eq.(2.6) we have £0 — (/3/2i^)6oc(0o)- From the 
second line of Eq.(2.15), ip(z) is independent of tpo and 



decays slowly as z for z ^ Iq. On the other hand, 
the local free energy density /loc decays as 



hoc/kBT,^2{P/2vYA,{z^l^) 



(2.17) 



In the right hand side. A,, is the universal number in 
Eq.(2.8) and hence the coefficient in front of (z -t- 
is universal. If we assume the linear form fs = ~hi^{Q), 
where hi is called the surface field. From the profile 
(2.15) hi and i/'o are related by hi = {(3 /v)C{tpo)4>o/ 

On the other hand, in the off-critical case ipoo 7^ 0, 
il>{z) approaches V'oo exponentially for large z as 



i^(z) - i/'cx) ~ iV'ool exp(-z/^oo)- 



(2.18) 



We introduce the correlation length = (C//")^/^ in 
the bulk region at -0 = ipoo, where /" = d'^f/d'^'4>. From 
Eq.(2.3) it is calculated as 



Coo = [C^/5{5+\)B^Y/^\i,oo\-'''^ 
= 6oc(V'oo)/V2'5((5 + l) 



(2.19) 



The second line is written in terms of ^loc in Eq.(2.6). In 
Fig.l, we show the scaled profile tl'{z)l\^oo\ vs z/^oo at 
-00 — 201-000 1- For small -000, 0(z) approaches \ipoo\ at 
z ~ Coo- For -000 > 0, the changeover from the algebraic 
decay to the exponential decay then takes place. For 
"000 < 0, ^{z) further changes from positive to negative 
on the scale of Coo- The length Coo is proportional to 
|'0oo|~''^'^, so it becomes longer with decreasing |'0oo|- 

In the off-critical semi-infinite case, the excess adsorp- 
tion Fsomi = /q dz(%p{z) — ipoo) is finite, where the sub- 
script ± represents the sign of ijjoo ■ Numerically we find 



scm± 



O 



(^AD/^o)V2fc 



(2.20) 



where k ^ i^/P - I, B+ = 0.310, S_ = 1.043, and iJjd is 
defined in Eq.(2.9). This formula includes the correction 
due to finite 0o and, as a result, it holds within 0.1% 
for |V'o/'0oo| > 20. We also recognize that Fsom- is three 
times larger than rsom+ for the same |'0oo|- The excess 
adsorption is larger for tpoo < than for tpoo > 0. In this 
relation and those to follow, we may push -00 to infinity 
to obtain the asymptotic relations near the criticality. 

So far, 00 is assumed to be small, so the transition 
length is given by £0 in Eq.(2.16). For large 0o of order 
unity, 0(z) should decay into small near-critical values if 
z exceeds a microscopic distance. 



C. Profiles between parallel plates 

We assume that a mixture at T — Tc is inserted be- 
tween parallel plates separated by D. The plate area S 
is assumed to be much larger than such that the edge 
effect is negligible. The fluid is in contact with a large 
reservoir containing the same binary mixture in equilib- 
rium. In the reservoir, the mean order parameter is "000 
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FIG. 1: (Color online) Normalized order parameter tp[z) /\'il)oo\ 
vs z/^ao at T = Tc in the semi-infinite case at i/jq = 20|'!/jcx) 
for positive and negative small ipoo- The decay is algebraic 
for small z as in Eq.(2.15) and is eventually exponential as in 
Eq.(2.18). For ipao < 0, the approach to tpoo is slow and the 
adsorption is large as in Eq.(2.20). 



and the chemical potential is /ioo in Eq.(2.10). Then we 
minimize the excess grand potential (per unit area), 



D 



dz[f loc - /oo - Moo (V" - i>oo)] 



(2.21) 



The fluid is in the region < z < D. At the walls z — 
and D, we assume the symmetric boundary conditions, 

^(0) = ^0, i^{D) = V'o, (2.22) 

where tpo > 0. Here Eq.(2.11) still holds, leading to 

C{ijj)\d^/dz\'^ = 2uj, (2.23) 

where cj is the excess grand potential density for a film, 

^ = /('/')- /m - Moo (V' - V'm) 

= + /,„ - /oo - Moo(^/'m - V'oo) (2.24) 

Hereafter, 4'm = 4'{D/2) is the order parameter at the 
midpoint z — D/2 and /,„ — /(V'm)- From Eq.(2.13) 
uj — LJs is simply a constant. We require dip/dz = at 
z = _D/2 or w = at "i/i = ipm because of the symmetric 
boundary conditions in Eq.(2.22). In the region < z < 
D/2, -ip = ^{z) is obtained from 



D/2~z^ r #VC(V')/2w(V'), 



(2.25) 



where a; = u;('0) is treated as a function of V'- As z — >■ 0, 
the plate separation distance D is expressed as 



D = / dV'\/2CMA^, 

'ipm 



(2.26) 



which determines ipm for each D, indicating the follow- 
ing, (i) In the integral of Eq.(2.26), we may push the 
upper bound ipo to infinity for large ipo since 1 / ^ 
^-(i+'5)/2 fQj. large ip- Thus V'm becomes independent of 
"i/jo as "i/jo OO- (ii) For a thick film with D > ^oo, V'm 
should approach ipooi where ^oo is the correlation length 
in Eq.(2.19). In the integrand of Eq.(2.26), we expand lu 
in powers of ip = ip — ipm as 



W = (/im - Moo)<(5 + /"(V'm)</'^/2 + 
= /"(V'm)[(V'm-V^oo)</^ + </^V2], 



(2.27) 



where /i,„ = fj,{ipm)- The second line follows for ipm — 
-0OO- Using ^oo = [C{ipoo)/f"{ipoo)]^^'^, we perform the 
integral in Eq.(2.26) as 2^oo ln[V'oo/ (V'm - "(Aoo)] for ipoo 7^ 
0. Thus, for D > ^00, 4'm approaches "000 as 



'ipni-lpoo ^ Ipoo exp(-i:»/25oo)- 



(2.28) 
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FIG. 2: (Color online) Normalized order parameter ip[z)/ipD 
vs z/D at r = Tc in a film in the half region Q < z < D /2 
for s — ipoo/tpD ~ 2,1,0,-1.11, and -1.5 from above with 
ipo ~ 20ipn in the symmetric boundary conditions, where ipn 
is defined by Eq.(2.9). The decay is algebraic for small z as 
in Eq.(2.15). 



In the limit ipo —> co, the profile ip{z) is scaled in terms 
of a scaling function 5'(w, s) as 



V'(z) 
Hereafter, we set 



lPD'^{z/D,lPao/'lpD)- 



S = 1poo/lpD, 
m = V'm/V'D = *(0.5, s). 



(2.29) 

(2.30) 
(2.31) 



which are the normalized order parameter in the reser- 
voir and that at the midpoint, respectively. In Fig. 2, for 
'0o/'0-D — 20, we plot ip{z)/ipD vs z/D for five values of 
s. These profiles represent '^{z/D,s) slightly away from 
the wall or for z ^ Iq. For s = 2 and —1.5, we find 
m = s (or 1pm — Ipoo)- However, the curves for s = and 
s — 1 are very close, while that for s = —1.11 tends to 
at the midpoint. 
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FIG. 3: (Color online) Left: m — tl^m/tpD vs s = tpao/tpD for a 
film at T = Tc in the limit tpo ^ oo, where i/'m is the midpoint 
value, ipaa is the reservoir value, and tpo (x D~^^" . Here m 
stays at 0.8 in the weak response region (—0.9 < s < 0.7), 
changes steeply in the catastrophic region ( — 1.2 < s < ~0.9), 
and tends to s in the strong reservoir region (|s| > 1.1). Right: 
Normalized adsorption T/DtpD vs s sd, T — Tc- It approaches 
2B+s~^ for s > 1 and 2B_|s|-''' for s < -1 (dotted lines) 
with k ^ v/P -1 from Eq.(2.20). The lines of T/Di/'d and 
2B-\s\^^ crosses at s = —0.90, where the Casimir amplitude 
A(s) is maximized from Eq.(2.52). 



D. Weak response and catastrophic behaviors 

In the left panel of Fig. 3, we show m vs s determined 
by Eq.(2.26) in the limit t/jq ^ oo- The slope of the curve 
dm/ds is written as 



dm 
ds 



&tpr, 

dtp, 



oo / D 



Xrj 



(2.32) 



which is equal to the ratio of the two susceptibili- 
ties Xm — kBTcdipm/ d^aa at the midpoint and Xoo = 
kBTcdtpao/dfioo in the reservoir—. Here, because of the 
large size of the critical exponent S = 4.815, there are dis- 
tinctly different three regions of s: (i) For —0.9 ^ s < 0.7, 
m is nearly a constant about 0.8 with dm/ds ~ 0.1. 
Here m little changes with a change in s. The response 
of t/j in the film to a change in ipoo is weak, (ii) For 
— 1.2 ^ s < —0.9, we have dm/ds ~ 5 ~ (5. In this catas- 
trophic region, m changes steeply between the base curve 
m — s and 0.8. The width of this region is of order 1/6 
and is narrow, (iii) In the regions s ^ 1 and s < —1.3, 
Eq.(2.28) gives 



exp[ 



(2.33) 



That is, iprn — "000 and the reservoir inflence is strong. 

The excess adsorption F = dz{ip — ipoo) in the film 
with respect to the reservoir is expressed as 



^0 



dV'(^-V'oo)v/2C(V)/w(^) 



(2.34) 



Again we may push the upper bound of the integral to 
infinity for large ■(/'o; then, T/DtpD becomes a universal 
function of s. At T = T^, Eq.(2.20) indicates that F ^■ 



2r,em± = 2DiI;dB±\s\^-''/^ for -D > Coo or for |s| > 1. It 
is convenient to introduce a scaling function F*(s) by 



F = 2F,em± + DtPoT* 



(2.35) 



In the right panel of Fig. 3, we show T / DipD vs s in the 
limit tpQ — >• oo. In the weak response region, —0.9 ^ s ^ 
0.7, we have T /DtpD ~ 0.5 — s. Thus T/DtpD increases 
with decreasing s, exhibiting a peak at the border of 
the weak response and catastrophic regions. In fact, its 
maximum is 2.29 at s = —0.94. In the strong reservoir 
region, it approaches 2B±\s\^'^ with k — v/jS — l (see the 
dotted hnes in the right panel of Fig. 3). 



E. Casimir term in the force density 

Using Eq.(2.23), the excess grand potential VI in 
Eq.(2.21) is expressed in terms of cj = w('0) as 



+2 / di>^/2CW)^- 



(2.36) 



In the right hand side, the first term arises due to the 
reservoir. In the second term, the gradient contribu- 
tion gives rise to the factor 2. Let us calculate the 
derivative dVt/dD at fixed and ^oo, treating ipm as 
a function of D. The derivative of the second term 
is D{poo — p^m)dipra/ dD from Eq.(2.26), where /x™ = 
/'(V'm)- We then find a simple expression, 



dD 



frn foo f-^oo (V^m V'oo ) ■ 



(2.37) 



The Casimir amplitude of the force density A 
{dft/dD)D'^ / ksTc is expressed as 



A 



D" 



ksTc 



fm foo /^oo 



(2.38) 



Note that the osmotic pressure H is the force density per 
unit area exerted by the fluid to the walls. See Appendix 
A for more discussions. Thus we fincti^ 

U^~—^^kBT,A/D''. (2.39) 
oU 

From the second line of Eq.(2.24) we also notice the rela- 
tion uj{ip) - ujs{ip) = -dn/dD. Equations (2.36)-(2.39) 
hold even for T ^Tc 'va our theory in the next section. 
At T = Tc, use of Eqs.(2.3), (2.8), and (2.9) gives 



A, 



A{s)=6\s\'+^ + \m\'+^ -{l + S)\s\'- 



(2.40) 



in terms of s in Eq.(2.30) and m in Eq.(2.31). We can see 
that A = A{s) is a universal function of s as V'o ~^ oo- 
In Eq.(2.26) we set ip — ipoq to obtain 



dqq-^''/^f'/^u;{q,s,m) = 1, 



(2.41) 
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where uj = uj/BqiP]^^ depends on q, s, and m as 



- 17711^+" - (1 + <5)|s|''-\s(q - to). (2.42) 

We seek to = m(s) for each s from Eq.(2.41) (as given in 
the left panel of Fig. 3). In particular, for \s\ ^ 1, we have 
"0™ — >■ ipoo and dfl/dD oc (i/'m — "ipoo)^ from Eq.(2.37). 
From Eq.(2.28) A{s) decays for \s\ > 1 as 



-4(s) 



exp 



hv/2^(J + l)|sr/'^]. (2.43) 



We calculate A{s) for two special cases, (i) First, for 
s = 0, the reservoir is at the criticality. Here, to > 0, so 
setting V = q/m in Eq.(2.41) gives to"/^ — /q with 



(2.44) 



where /q = 0.719 for = 3. Since dVl/dD = f{ipm), we 
obtain the critical-point value ^cri = -^(0) in the form. 



A ■ — A 

-^cri — c ; 



(2.45) 



where Ac is given by Eq.(2.8) and will be estimated be- 
low Eq.(3.16). Essentially the same calculation was origi- 
nally due to Borjan and Upton^ (ii) Second, we assume 
'4'm = 0, which is attained for ipoo < or for s < 0. See 
the corresponding curve of s = —1.11 in Fig. 2. From 
Eq.(2.41) we obtain |s|''/'^ = Jo with 



Jo 



dvv-^T'''^^ /^v^+^ + {l + 5)v, (2.46) 



where Jq = 1.195 at d = 3 numerically, leading to s = 
-J^'" = -1.1170. From Eqs.(2.40) and (2.45) we find 

^(-1.1170) = JiAcS = AcMo/Iof6, (2.47) 

Thus yl(-1.1170)/Ari = 24.45 for d = 3. 

In Fig. 4, we display .4(s)/^cri for d = 3 calculated 
from Eqs.(2.40) and (2.41). Its maximum is 24.45 at 
s = —1.1170, where ipm — 0. To be precise, the curve 
exhibits a small cusp due to the weak -^-dependence of 
C{tp)- Similar enhancement of A was found at off-critical 
compositions by Maciolek et for two-diemnsional 
Ising films and by Schlesener et al^ in the mean-field 
theory at T = Tc. In the next section, the origin of this 
peak will be ascribed to the fact that the peak point on 
the line T = Tc is very close to a capillary-condensation 
critical point in the region T < Tc (see Figs. 9 and 12). 

Mathematically, the peak of A{s) at T = Tc stems 
from the presence of the weak response and catastrophic 
regions, for which see the explanation of the right panel 
of Fig. 3. Here we calculate dA{s)/ds from Eq.(2.40) as 



dA 



as- 



-1 ds 





TO 


(5-1 

TO 




( 


S 


T " 






- 28Ac - 


- 40 



dm 
ds 



— TO + s. 



(2.48) 



region, we may neglect the first term in the right hand 
side to obtain dA/ds = —a\s\^~^{m — s) < 0, which is 
nearly zero for s > —0.5 and grows abruptly for s < —0.5. 
In the catastrophic region, dm/ds is of order S such that 
dA/ds changes its sign, leading to a maximum of A. 



F. Casimir term in the grand potential 

For binary mixtures, the de Gennes-Fisher scaling 
formJ'^ for the grand potential reads 



n = n^- knTD-'^'^-^^A, 



(2.49) 



where Qcx, = lim£)_>oof^ is the large-separation limit. 
The A is a function of s in Eq.(2.30) (and a scaled 
reduced temperature in the next section). This form 
may be inferred from the boundary behavior of /loc in 
Eq.(2.17). Next, we differentiate Eq.(2.36) with respect 
to /Zoo (or V'oo) at fixed D and ipo- Following the pro- 
cedure used in deriving Eq.(2.37), we obtain the Gibbs 
adsorption formulapii. 



an 

9/ioo 



Xc 



on 



ksTc dtpo 



-r. 



(2.50) 



where F is the excess adsorption in Eq.(2.34) and Xoo — 
kBTc/f'ii^ac)- If the Fisher-de Gennes form (2.49) is 
substituted into Eq.(2.50), the above relation yields 



a7 



— ^ooF* 



(2.51) 



where F* is defined by Eq.(2.35) and Roo is the follow- 
ing dimensionless combination representing the scaled in- 
verse susceptibility in the reservoir. 



Roc = D'^^l/Xoo. 



(2.52) 



At T Tc, we obtain R^o = Ac5{5 + l)|s|''"\ which has 
already appeared in Eq.(2.48) as olsl"*^^. In our theory 
in the next section, Eqs. (2.49)-(2.52) will remain valid 
even for T 7^ T^. 

At T = Tc, direct differentiation of 17 in Eq.(2.49) with 
respect to D also yields a relation between A{s) and A(s), 



^(,s) = (d-l)A(.)-^^A(.). 



(2.53) 

Elimination of dA/ds from Eqs. (2. 51) and (2.53) yields 

(2.54) 



A{s) = 4^ + 4^s\srhr*is). 



d-1 d-1 

Thus, as well as A, A is proportional to the universal 
number Ac. The critical-point values of the amplitudes, 
written as Acri — A(0) and Aa-i — A{0), are related as 



Ac 



Acri/ {d- I). 



(2.55) 



In Fig. 4, we also plot A(s)/Acri for d = 3 numerically 
calculated from Eq.(2.54). Its peak height is 13.36 at 
s = —0.90, which is about half of the height oi A{s)/Acii- 
The differential equation (2.53) is excellently satisfied by 
A{s) and A(s) in Fig. 4. From Eq.(2.51) the amplitude 
A(s) is maximized at a point where F* — 0. Indeed, in 
the right panel of Fig. 3, the line of T/DipD vs s and that 
of 2B_|s|''/^^^ vs s cross at s = -0.90, where A(s) is 
maximum in Fig. 4. 
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FIG. 4: (Color online) Casimir amplitude ratios A{s)/Acii 
for the force density and A(s)/Acri for the grand potential vs 
s = tpoo/tpD at T = Tc in the limit ipo oo. These quantities 
are maximized at s ~ — 1 upon the changeover between the 
weak response and catastrophic regions (see the text). 



III. CRITICAL BEHAVIOR FOR T ^ n 

In this section, we now set up a local functional the- 
ory including the gradient free energy for nonvanishing 
reduced temperature, 

T = T/Tc - 1. (3.1) 

For binary mixtures, we suppose the upper critical solu- 
tion temperature (UCST), while r should be defined as 
T = 1 — T/Tc for the lower critical solution temperature 
(LCST). Our model is similar to the linear parametric 
model by Schofield et a^42r— and the local functional 
model by Upton et aL— '2^'^'^. The latter is composed 
of the free energy of the linear parametric model and the 
gradient free energy. We use a simpler free energy ac- 
counting for the renormalization effect and the gradient 
free energy. Further, we define the free energy within the 
coexistence curve and make it satisfy the two-scale-factor 
universalit y '^^i'^^ . Appendix B will give the relationship 
between our model and the linear parametric model. 

The Casimir amplitudes A in Eq.(2.38) and A in 
Eq.(2.49) depend on s = ipoo/i^D in Eq.(2.30) and the 
scaled reduced temperature, 

t - T(i?/eo)^/^ (3.2) 

where is a microscopic length (^ 3A) in the correlation 
length ^ = ^qt^"^ for r > at the critical composition. 
For example, ^/D is 0.24 for t = 10 and 0.056 for t = 100. 

A. Model outside the coexistence curve 

For T < 0, two phases can coexist in the bulk. The 
coexistence curve is written as '0 = iV'cx with 

V-cx = b,^\Tf, (3.3) 



where &cx is a constant. In this subsection, we present a 
local free energy density f\oc applicable outside the coex- 
istence curve (IfAl > V'cx if T < 0). In the next subsection, 
we present its form within the coexistence curve. 

As a generalization of the model in Eqs.(2.1)-(2.4) for 
T = 0, we again use the usual form, 

fioc = f+\kBT,C\Wij\\ (3.4) 

The free energy density / — f{ip,T) is an even function 
of ijj expressed as 

/ = fcsr,Qr^2 + iu^4y (3.5) 

Here r, u, and C are renormalized coefficients depending 
on ^ and r. To account for the renormalization effect, 
we introduce a distance from the critical point w. In the 
critical region (w <C 1), these coefficients depend on w as 

C = Ciw-'T', (3.6) 
r = C^^^w-'-'t, (3.7) 
u = Cl'^i*eo"^w;(^~2^)^ (3.8) 

where Ci and u* are constants. For r > and i/j = 0, 
we set w = T to obtain — C jr. For t = 0, we should 
have up'^ cx ■0^ . We thus determine w as a function of r 
and i]?' from 

u; = r + C2W^"2^V^ (3.9) 

where Ci is a constant. Recall the mean-field expression 
— r + iwp'^ for the susceptibility Xj which holds even 
for r < 0. If we use Eqs.(3.6)-(3.8), we have r + 3utp^ = 
Ci£,o'^w''~^[t + 3Ciu*^o"%i-2^'0^] in our renormalized 
theory. Thus we may require 

r + 3wjj^^Ci^Q^w'>, (3.10) 

where we relate C2 in Eq.(3.9) to Ci as 

C2 = 3u*CiCo~'- (3-11) 

See Fig. 5 for / and w as functions of ip for fixed posi- 
tive T. In the simplest case r > and -0 = 0, we have 
w = T and w'^ — ^o/C so that the concentration suscep- 
tibility x{tj ■0) — ksTcidtl) / d^)r grows strongly with the 
exponent 7 as 

x{r.Q)^r^C^'eoT-\ (3.12) 

In the renormalization group theory*^, Eqs.(3.6)-(3.9) 
follow if iQ^w" is the lower cut-off wave number of the 
renormalization. The coupling constant u* should be a 
universal fixed-point value. Its e expansion reads 

Kdu =e/9 + ---, (3.13) 

where K^i is the surface area of a unit sphere in d dimen- 
sions divided by {2t:Y (so = l/87r2 and = l/27r2). 
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Retaining the small critical exponent r] oa and using 
the relations among the critical exponents, we may cor- 
rectly describe the asymptotic scaling behavior, though 
the critical amplitude ratios are approximate. In addi- 
tion, note that a constant term independent of has been 
omitted in the singular free energy density in Eq.(3.5), 
which yields the singular specific hea.t^^. 

For r = 0, the Fisher-Au Yang model in Eq.(2.1) folows 
with Co and Bq given by 



(e/2-r,)z.//3 



(3.14) 
(3.15) 



Thus we have BokBTjC'^ = u'C"^ ^"^^ I and Bq/Co = 
M*CiC2^^~^fcr74 = C2/'^/12^2_ rpj^g universal number 
Ac in Eq.(2.8) is calculated as 



A, 



(3.16) 



As a rough estimate, we set u* = l/9i^3 = 2Tr^/9 for d — 
3fromEq.(3.13). This leads to = {2^lV&)lu* = 1.49. 
This value yields Acri = Ari/2 = 0.279 from Eq.(2.45). 
For d = 3, KrecbSl estimated Acri to be 0.326 by a 
field-theoretical methods and 0.345 by a Monte Carlo 
method, Borjan and Upton= obtained Acri = 0.428 by 
the local functional theory, and Vasilyev et alM- found 
Acri = 0.442 by a Monte Carlo method. 

To make the following expressions for / and /i simpler, 
we introduce the dimensionless ratio, 



S — t/w, 
in terms of which -0^ is expressed as 

^2 = (l-5)u-2^/C2. 

Then / and ^ — {df /dip)T are expressed as 
/ 



w 



2-a 



ksTc 



{1 + ?>S){1- S), 



(3.17) 



(3.18) 



(3.19) 



^-" + <^7")^ + ^"f.^C.0. (3.20) 



In Eq.(3.20), we have used the relation {dw/d'ip)r — 
2C2W^-'^^^l[2fi 2/3)5']. However, the susceptibil- 

ity x{t^^) = ksTc/idfi/d^) is somewhat complicated^. 
It can be simply calculated for = as in Eq.(3.12) for 
T > and as in Eq.(3.24) below for r < 0. 

We now seek the coexistence curve (3.3) by setting 
11 = with T < 0. From Eq.(3.20) it follows the quadratic 
equation 2 — a + 4(1 — a)S + 5aS^ = of S*, which is 
solved to give 5 = — 1 /cr or 

w = a\T\. (3.21) 

The coefficient a is expressed in terms of a as 

cr = 2 - 9a/[2 + - 18a + 9a^], (3.22) 



which is equal to 1.714 for d = 3. If we substitute 
Eq.(3.21) into Eq.(3.9), the coefficient 6cx in Eq.(3.2) is 
calculated as 



(3.23) 



leading to C2 = {l+<7)a^P~^ /bl^. Since 6cx is experimen- 
tally measurable, there remains no arbitrary parameter 
with Eqs.(3.11) and (3.23). The susceptibility x on the 
coexistence curve is expressed as 



xi-h\,A 



3(7 



2C(l + a 

In the denominator of this relation, we introduce 
2(1 — a)cr — ba 



c 



(2/3(7 - 1 + 2/3)2 ' 



(3.24) 



(3.25) 



which is equal to 4.28 for d = 3. The ratio of the suscep- 
tibility for r > and -0 = and that on the coexistence 
curve at the same Irl is written as 



X(|r|,0) 2 -^^W7-i 

x(-rl,V'oo) 3 



(3.26) 



which is 8.82 for d = 3. Note that the e expansion 
gives i?^ = (2 -I- e)2'/^ + • • • and its reliable estimate 
is 4 943, 52, 53 ^ g^jgQ ^Yite the correlation length on the 
coexistence curve as ^ = (C//")^^^ — '^o!'''! where 
^0 is another microscopic length. The ratio of the two 
microscopic lengths ^0 and is written as 



(3.27) 



which gives Co/^o ~ 2. 99 for 0? = 3. Note that the e 
expansion result is £,o/£,q = 2''(1 H- 5e/24 -I- • • •) and its 
reliable estimate is 1.9 for d = 343,52,53 ^ recognize 
that the correlation length ^ and the susceptibility x on 
the coexistence curve are considerably underestimated in 
our theory (mainly due to a factor Z{9) in Eq.(B13) in 
Appendix B). 

Finally, let us consider the characteristic order param- 
eter of a film V'c defined in Eq.(2.9). From the sentence 
below Eq.(3.15) it is written as 



4,D = C^"\V2iS,o/Uf/'' 
= 1.47fecx(eo/^)^/" 



(3.28) 



where 6cx is given by Eq.(3.23). Equation (3.18) gives an 
expression for ip valid outside the coexistence curve, 

W\/i^D = 24-'3/2-(i _ 5)i/2|5p/J|i|/3_ (3_29) 
On the coexistence curve, this expression becomes 

V-cx/^u = 24-^/2^(1 + a-i)i/V^|i|'3, (3.30) 
which is equal to 0.66|ip in our theory. 
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l/'/V'cx V'/'Acx 



FIG. 5: (Color online) Normalized singular free energy density 
f{ip)/fcri (left) and normalized distance from the criticality 
w/\t\ = 1/|5| (right) vsV/iAcx, where /cri = kBn\Tf''°'/u*^i 
and Tpcx = bcxl'rl'' (both for r > and r < 0). The upper 
(lower) curves correspond to those for r > (r < 0). For 
T < 0, w is independent of ip within the coexistence curve 
liAI/iAcx < 1 in our model. 



are continuous only up to ^ = 2 on the coexistence curve. 

Within the coexistence curve, we assume a 'ip^ theory 
including the gradient free energy. Since fi = df/dip 
vanishes on the coexistence curve, / is of the form, 

/ = /cx + IkBTMi'^ - C)', (3.31) 

where /cx is the value of / on the coexistence curve. From 
Eq.(3.19) it is written as 

/cx = kBT,{a - 5)(1 + (7)|r|2-V36u*eoa". (3.32) 

We determine the coefficient i?cx requiring the continuity 
of the second derivative /" = fcsTc/x- Then we obtain 
2_BcxV'cx = where x is given by Eq.(3.24). Some 

calculations give a simple result, 

Bex = C"cx, (3.33) 



4 





/^oo = 




i;^^ t = -15 ('^oo = ■'/'cx) 




■ --— (V'oo = V'cx) 
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FIG. 6: (Color online) Normalized order parameter ^p{z)/tpD 
vs z/D for = in a film in the half region < z < D/2 
for t = 30, 15, 5.51, -4, and -15 from below with ipo ~ 20i/'i5 in 
the symmetric boundary conditions. Here ipoo = for t > 
and tpao = V'cx for t < 0. The decay is algebraic for small z 
as in Eq.(2.15). 



B. Model including the coexistence curve interior 

For T < 0, we need to define a local free energy density 
/ inside the coexistence curve IV"! < V'cx, where the fluid 
is metastable or unstable in the bulk. Notice that ipiz) 
changes from — V'cx to V'cx in the interface region in two- 
phase coexistence. Thus, to calculate the surface tension 
in the Ginzburg-Landau scheme, we need /(V") in the re- 
gion IV' I < V'cx for T < 0. Since the linear parametric 
model is not well defined within the coexistence curve. 
Fisher et alM proposed its generalized form applicable 
even within the coexistence curve to obtain analytically 
continued van der Waals loops. We propose a much sim- 
pler model, though the derivatives d^f/dip^ in our model 



where ^ is given by Eq.(3.25) and Ucx is the value of u in 
Eq.(3.8) on the coexistence curve written as 

Mcx = C2u*e-^|aT|(^-2^)^ (3.34) 

It follows the relation /cx — ksTcUcx'^'txi'^ ~ 5)/4(l + a). 
The coefficient of the gradient term is replaced by its 
value Ccx on the coexistence curve: 

Ccx = Cl|ar|-''^ (3.35) 

Then the susceptibility % and the correlation length ^ are 
continuous across the coexistence curve. In this model, 
the renormalization effect inside the coexistence curve is 
the same as that on the coexistence curve at the same r. 
See Fig. 5 for / and w vs ip for fixed negative r. 

The surface tension (t^ is given by the classical formula 
(Ts = 2/cbTcCcxV'cx/3'C in the (/)^ theory^. It is propor- 
tional to ^-''+1 with ^ — Cokl Some calculations give 
the universal number. 

As = <JsC''''/kBT, 

= [2C(1 + (T"^)/3]^/V3C"*- (3.36) 

Using Eq.(3.25), we have Ag = 0.165/u*. Further, if we 
set u* = 27r2/9, we obtain Ag — 0.075. Note that A^- is 
known to be about 0.09 for Ising-like systems^. 



C. Casimir amplitudes for T Tc 

From Eq.(2.38) we can readily calculate A{t, s) numer- 
ically as a function of s and t. However, to calculate 
A(i, s), we cannot use Eq.(2.54) for t ^ and need to 
devise another expression. To this end, we write the sec- 
ond term in the right hand of Eq.(2.36) as 

J = n- D— = 2 d^V2C(V')w(V'), (3.37) 
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FIG. 7: (Color online) Left: Normalized midpoint order pa- 
rameter m{t) = tpm{t)/ipD vs t = T^D / S^oY^" for a film with 
thickness D, where s = for r > and s — ■(/'cx/'f/'cx for r < 
(Moo = 0). The coexistence curve tpcx{t)/tpD vs t in Eq.(3.3) 
is also shown in the region t < 0. Right: Normalized excess 
adsorption T/DiI^d vs t for fioo ~ 0, whose maximum is 1.51 
att^O. It tends to 2g±\tf~'" for t>W and for t < -1 with 
g+ = 1.245 and g- = 0.471, respectively (see Eq.(3.45)). 



where is defined by Eq.(2.24) with / being given 

by Eqs.(3.5) and (3.31). In the limit D — ^ oo, we have 
J — >■ iloo, where fioo is the large-separation limit: 



-00 



(3.38) 



Here, the lower bound of the integration is i/'oo and a;s(V') 
is the grand potential density for the semi-infinite case 
in the form of Eq.(2.13). Dividing the integration region 



and [ipao , "0™] ) we find 



in Eq.(3.38) into [ipm,ipo] 

-2 / d'^p^J2C{^l))uJs{1l')■ (3.39) 

In the first integral, we may push the upper bound ■i/'o to 
infinity, since the integrand tends to zero rapidly for large 
tp from -y/cJ — ^JUl — —[d^l/dD)/[^Ju + ^JUl] (see the 
sentence below Eq.(2.39)). From Eqs.(2.37) and (2.49) 
we obtain 



D 



■(J-f^oc). 



(3.40) 



With this expression, we can calculate A(s,t) numeri- 
cally, We confirm that it is a function of s and t only. 

In our theory, Eqs.(2.36)-(2.39) and Eqs.(2.49)-(2.51) 
remain valid even for t ^ 0. We now treat to. A, A, and 
r* as functions of s and t, so d{- ■ ■)/ds = {d{- ■ ■)/ds)t 
and d{- ■ ■)/dt = {d{- ■ ■)/dt)s. From Eq.(2.38) we obtain 
the generalized form of Eq.(2.48) as 



1 dA 

Roo ds 



Mo 



-Xc 



dm 
'97 



(to — s). 



(3.41) 



Here, Roo is a function of s and t defined by Eq.(2.52) (see 
Fig. 16 for its overall behavior in the s-t plane). Although 
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FIG. 8: (Color online) yl(i)/Ari vs t = r(D/^o)^'''' along 
the critical path /ioo = from our theory (red bold line) and 
from Borjan and Upton's theorj^ (blue broken line), which 
are in good agreement with the Monte Carlo data (-f)'^^. Also 
plotted is A(t)/Acri vs t from our theory on the critical path. 



T>T, 




FIG. 9: (Color online) Normalized midpoint order param- 
eter m{s,t) = ^m/tpD (top left), normalized excess ad- 
sorption T/tpoD (top right), A{s,t) / Acri (bottom left), and 
A(s,t)/Acri (bottom right) for f > in the s-t plane. The 
latter three quantities are peaked for s ~ — 1. The amplitudes 
A{s,t) and A{s,t) are very small for t ^ 1 and \s\ 2> 1. 



redundant, we again write Eq.(2.51) as 



Roo ds 



(3.42) 



which follows from the Gibbs adsorption relation. Here, 
r* is defined by Eq.(2.35), where the excess adsorption in 
the semi-infinite case Fgeni-i- should be calculated for each 
given ipoo and r. In addition, differentiation of Eq.(2.49) 
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with respect to D yields the generahzation of Eq.(2.53): 

^=d-lA-^- — . 3.43 

V OS V at 



D. Results for T > 

We first give some analysis along the critical path 
Moo — 0, where V'oo = for i > and il^oo — V'cx > for 
t < 0. Numerical and experimental studies on the crit- 
ical adsorption and the Casimir amplitudes have mostly 
been along this path in the literatur o^'^i^^i^''~^^i'^^i^^i^^ . 

Since Eqs.(2.23)-(2.26) stiU hold, ipiz) may be calcu- 
lated in the same manner as in Borjan and Upton's paper 
on the critical adsorption^. In Fig. 6, we show ip{z)/'4iD 
for various t at ipo/ipD = 20. For a film with large posi- 
tive t, Eq.(2.26) gives m = 'ipm/4'D in the form, 



E. 



eM-D/20 exp(-r/2). 



(3.44) 



For t < -1, Eq.(2.28) gives m - s - sexp(-D/20 = 
sexp(-^o|ir/2Co) < 1. In Fig.7, we plot m(t) = 
tpm{t)/ipD and r{t)/DipD vs t. We can see that m(t) 
decays as in Eq.(3.44) for r ^ 1 and tends to ipcx{t)/4'D 
for r <C — 1. As discussed around Eq.(2.35), r{t) tends 
to 2rscm± (t) , where Fsomi is the excess adsorption in the 
semi-infinite case. In our theory Fgemi behaves on the 
critical path as 



13-,. 



g±\tf--D^l, 



(3.45) 



where g+ = OMg'+ = 1.245 for i > and = OMg'_ = 
0.471 for i < 0. We obtain the ratio g+/g- = 2.64, while 
it was estimated to be 2.28 by Floter and Dietrich^. The 
right panel of Fig.7 shows that F approaches the limit 
2rsom+ for t > 10 and 2rsom- for i < -1. For < < 0, F* 
in Eq.(2.35) is very small. For example, it is —0.0675 at 
t = -1 and -0.0122 at t = -2.2 . 

Figure 8 displays the normalized amplitudes A{t) /AcH 
and A(t)/Acri in our theory. We also plot A{t)/AcTi from 
the Monte Carlo calculation by Vasilyev et al^ and from 
the local functional theory by Borjan and Upton^^. Re- 
markably, the two theoretical curves of ^(t)/^cri excel- 
lently agrees with the Monte Carlo data in the region 
t > 0. We should note that the free energy density in our 
theory and that in the linear parametric model — used 
by Borjan and Upton^S are essentially the same on the 
critical path in the region i > 0, as will be shown in Ap- 
pendix B. In our theory, A/Acri and A/Acri behave sim- 
ilarly. The maximum of the former is 2.173 at t = 5.51, 
while that of the latter is 1.872 at t = 2.30. 

In the literature, however, there has been no calcula- 
tion of the Casimir amplitudes in the s-t plane for s ^ 
accounting for the renormalization effect. In Fig. 9, we 
display m(s, t), F(s, t)/^DD, A{s, t), and A(s, t) for i > 
in the s-t plane. We can see that A{s,t) and A(s,i) are 
both peaked at s ^ — I and behave similarly. 



Phase behavior of capillary condensation and 
enhancement of the Casimir amplitudes 
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FIG. 10: (Color online) Phase diagram of a near-critical fluid 
in a film for large adsorption in the s-t plane, where s = 
V'oo/'/'d(oc i/Joo-D''/'^) and t — t{D /^^qY^" . Two-phase region 
is written in the right (in gray). There appears a first-order 
phase transition line (red) of capillary condensation with a 
critical point at {s,t) — ( — 1.27,-3.14). Plotted also are a 
line of maximum of A{s,t) with {dA/ds)t — and a line of 
maximum of A{s,t) (broken lines). The former approaches 
the coexistence curve for t < —6 and the latter is very close 
to the capillary condensation line. On paths (a)-(g), m, F, A, 
and A are shown in Figs. 11 and 12. 



We now discuss the phase behavior in the region of 
t < and tpoo < 0. In Fig. 10, we show a first-order phase 
transition line in the region t < and s < outside the 
coexistence curve. As in Fig.ll, the discontinuities of the 
physical quantities across this line decrease with increas- 
ing t and vanish at a critical point {t,^oc) = {t^'^,'^T)- 
In agreement with the scaling theory^, we obtain 



-3.14(eo/i?)^/^ 
-1.27V'z5 = -1.876c 



(3.46) 

.{^o/Df'\ (3.47) 



We use the second line of Eq.(3.28) in Eq.(3.47). For r < 
T^^, let us write the transition point as ^oo = '0ca(T) or as 
s — Sca.{t) = iJjca.{T)/ipD- This line divides the condensed 
phase with m > in the range Sca < s < —ipcK/'tpD and 
the noncondensed phase with m < in the range s < Sca- 
In Fig.ll, the isothermal curves of m = il}{D /2) /'ipn and 
T/DtpD are shown as functions of s for t < 0, where they 
are continuous for t > tj:^ = —3.14 and discontinuous 
for t < f^^. It is the capillary condensation line for the 
gas-liquid transition^'^ and is also the two-dimensional 
transition line for Ising-like films^"— i^. 

In Fig. 10, we also plot two dotted lines. On one line, 
A{s, t) takes a maximum as a function of s at fixed t. For 
t < f^^, it is slightly separated from the transition line 
for t > —4.0 but coincides with the transition line for 
t < —4.0. On the other line close to the bulk coexistence 
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FIG. 11: (Color online) Normalized midpoint order param- 
eter m — ■(/;(_D/2)/?/>_D (top) and normalized excess adsorp- 
tion V/xpoD (bottom) vs s in a film with thickness D for 
t = -0.3, -2, -3.1, -4, -6, -8, and -10 from the right. See 
Eqs.(2.34) and (2.35) for V. For t < —3.1, these quantities are 
discontinuous across the capillary condensation line in Fig. 10. 



FIG. 12: (Color online) Casimir amplitude ratios 
A{s,t)/Acii (top) and A(s,f)/Acri (bottom) vs s for 
t = -0.3,-2,-3.1,-4,-6,-8, and -10 from the right, in 
the range s < —ipcx/ipD = — O.GGltj*^. Curves (c) corresponds 
to the capillary-condensation critical point. Those (d)-(g) ex- 
hibit a first-order phase transition, where A is discontinuous 
but A is continuous. 



curve, A(s, t) takes a maximum as a function of s at fixed 
t and we have {dA/ds)t = and T* = from Eq.(3.42). 

In Fig. 12, we plot the Casimir amplitudes A{s,t) and 
A(s,t) as functions of s for t < 0. For t < t'^'^, they grow 
very strongly in the condensed phase. As a marked fea- 
ture, A and A behave very differently for t < 0, though 
they behave similarly for t > in Fig. 9. These results 
are consistent with their derivatives with respect to s in 
Eqs.(3.41) and (3.42). Previously, enhancement of A was 
found close to the transition line in the condensed phase 
by Maciolek et alJ^ in two-dimensional Ising films and 
by Schlesener et al^ in the mean-field theory. 

In Fig. 11, the slope {dm/ds)t or the susceptibility 
dm/ds in Eq.(2.32) diverges as t is decreased to f^^. Thus, 
in Fig. 13, we plot {dm/ds)t vs s for t — —3.1 in the left 
panel and the curve of to vs s at i = —3.1 and -2 in the 
right panel. The curve oi t — —3.1 can well fitted to the 
following mean-field form, 

s-s^ = ^ca(i-<r)("i-<") + Sca(m-m^/)3, (3.48) 




-1.27 -1.265 -1.26 -1.255 -1.25 0.01 

s 



0.1 

\rn — nf^ 



FIG. 13: (Color online) Left: Susceptibility {dm/ds)t = 
{dtpm/ dtp oo)t at t = —3.1 slightly above the capillary- 
condensation critical point {f^^ = —3.14). Right: Relations 
between deviations s — s^^ and m — m'^^ near the capillary- 
condensation critical point a,t t — —3.1 and away from it at 
t — —2 on a logarithmic scale. The former can be fitted to 
the mean-field form (3.48). 



where (t^^, s^^, to^^) = (-3.14,-1.27,-0.172), 
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0.0425, and Bca = 0.11. This mean-field behavior near 
the capiUary-condensation critical point arises because 
the long wavelength fluctuations of tp inhomogeneous 
in the lateral plane have been neglected. The curve of 
t = -2is not weh fitted to Eq.(3.48) for \m-ml''\ > 0.3. 



F. Determination of the capillary condensation line 

In Fig. 14, we show the profile ip{z) at t = —4 for four 
values of s given by (A) -1.10, (B) -1.29, (C) -1.33, and 
(D) -1.50, where tpo = 20-0^- The two lines ip = ±V'cx are 
also shown, between which {ipl < Tpcx and the free energy 
density is given by the mean-field form (3.31). Here, be- 
tween (B) and (C), there is a first-order phase transition 
at s = —1.31, where the normalized midpoint value m 
changes discontinuously between 0.407 and —0.793. 







t = -4 
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z/D 

FIG. 14: (Color online) Normalized order parameter xp{z)/tpD 
vs z/D in a film in the half region < z < D/2 at t = —4, 
with ■j/jQ — 20iI;d in the symmetric boundary conditions. From 
above, (A) s = —1.10 (at maximum of A), (B) -1.29 (near 
maximum of A), (C) —1.33, and (D) —1.50. A first-order 
phase transition occurs between (B) and (C). Lines of ipao = 
ii/jcx are written, where tpc^/i'D = 1.04. 



In Fig. 15, A{s)/Aci-i, A(s)/Acri, and F* are displayed 
as functions of s at i = —4 (left) and -10 (right), where 
ip{z) for the points (A), (B), (C), and (D) at t = -4 can 
be seen in Fig. 14. In numerical analysis, we obtained two 
branches of the profiles giving rise to hysteretic behavior 
at the transition. From Eq.(2.49) A should be maxi- 
mized in equilibrium. Thus the equilibrium (metastable) 
branch should be the one with a larger (smaller) A. A 
first-order phase transition occurs at a point where the 
two curves of A{t) cross. We confirm the following, (i) 
In accord with Eq.(3.42), the maximum point of A co- 
incides with the vanishing point of F*, as can be known 
from comparison of the middle and bottom panels, (ii) 
As s -> — ?/'cx/'/'r>, A becomes very small on the top 
plates, since it approaches the value on the critical pass 
Moo = with t < displayed in Fig. 8 (see the sentence 
below Eq.(3.44)). (iii) As s is decreased from —ipc^/ipD 
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FIG. 15: (Color online) A{s,t)/Acri (top), A(s,t)/Acri (mid- 
dle), and r*{s,t) (bottom) for t = -4 (left) and t = -10 
(right). For t = —4, points (A), (B), (C), and (D) correspond 
to the curves in Fig. 14. There are equilibrium and metastable 
branches near the transition. Amplitude A is maximized at 
(A) (left) and at (a) (right), where F* = from Eq.(3.42). 
The transition occurs between (b) and (b'), where A is con- 
tinuous. 




FIG. 16: (Color online) Scaled inverse susceptibility for a film 
Roo in Eqs.(2.52) and (3.42) in the s-t plane, which is of order 

40151*"^ for \s\ > 1. 



14 



to the transition value Sca in the condensed phase, A 
grows with a steep negative slope. To understand this 
behavior, we compare the first and second terms in the 
right hand side of Eq.(3.41). We notice that the ratio 
of the first term to the second is very small in this s 
range. It is between [0.02,0.05] for t = —4 and between 
[0.001, 0.005] for t = —10. In the condensed phase in the 
range Sca < s < —ipcx/ipD, we thus find 

BA 

— ^-i?^(TO-s). (3.49) 
OS 

While m — s remains of order V'cx/V'-D ~ j^l'^i the co- 
efficient i?oo is very large outside the coexistence curve. 
Figure 16 displays the overall behavior of Roa in the s-t 
plane, where i?oo — 6.6|i|'*' on the coexistence curve and 
i?oo = AcS{l + 5)\s\^-^ - 40|s|*-i for |s| > 1. 

We also comment on the validity of the general relation 
(3.43) in our numerical analysis. For example, at t = —4 
in Fig. 15, let us consider the two points (A) s = —1.1 
and (B) s = -1.29. At point (A), we have A = 5.54, 
A = 9.02, dA/ds = 3.13, and dA/dt = -2.21. The 
three terms in the right hand side of Eq.(3.43) are thus 
18.04, 1.78, and —14.03 in this order and indeed their 
sum gives A. At point (B), we have A = 23.98, A = 1.98, 
dA/ds = 72.4, and dA/dt = -3.85, so the three terms 
in the right hand side of Eq.(3.43) are 3.96, —24 A, and 
48.2 in this order, whose sum indeed yields A. 

IV. SUMMARY AND REMARKS 

We have calculated the order parameter profiles and 
the Casimir amplitudes for a film of near-critical fluids. 
Our results are also applicable to one-component fluids 
near the gas-liquid critical point where the walls favor 
either of gas or liquid. We summarize our main results. 

(i) In Sec. II, we have used the singular free energy by 
Fisher and Au Yang aX T — Tc. Using this model, we 
have defined the two Casimir amplitudes, A for the force 
density and A for the grand potential, as functions of 
the scaled order parameter of the reservoir s = ipoo/'^D 
in Eq.(3.31), They are sharply peaked at s ~ — 1 and the 
peak heights are much larger than their critical-point 
values Acri and Acri as in Fig. 4. These singular behav- 
iors have been analyzed analytically. This off-critical 
behavior may also be interpreted as pretransitional 
enhancement, because the region of s ^ —1 and i = is 
close to the capillary-condensation critical point. 

(ii) In Sec. Ill, we have constructed a free energy with 
the gradient contribution for T ^ Tc including the renor- 
malization effects, with which we may readily calculate 
the physical quantities. The Casimir amplitudes A and 
A are much amplified for s ^ — 1 as shown in Fig. 9 for 
T > Tc and in Fig. 12 for T < Tc. Their maxima are 
larger than their critical-point values by 10-100 times. 
We have then found a first-order phase transition line 
of capillary condensation for negative t slightly outside 
the bulk coexistence curve, where the profile of ip and 



A are discontinuous but A is continuous. This line ends 
at a critical point given by {s,t) = (-1.27,-3.14). The 
amplitude A exhibits a maximum close to this line, while 
the amplitude A close to the bulk coexistence curve. 

We make some remarks. 

1) Even at the mean-field level, it follows the power-law 
form of the interaction free energy AF for ^ > D > Eq: 

AF = -keTcD-'^Ao. (4.1) 

Here, setting ^ = 3 in Eq.(2.3) and rj = in Eq.(2.4), 
we obtain Aq = AI^C^ /iB^ksTc at the criticality under 
strong adsorption, where /q is given by Eq.(2.45). The 
adsorption-induced interaction is already present in the 
mean-field theorjs^ii^ and its form becomes universal 
near the criticalitySS. Enhancement of AF near the 
capillary condensation transition is rather obvious in 
view of the fact that it occurs even in the mean-field 
theory^!. Note that AF cx h\ for weak adsorption 
(£o > D) away from the criticality, where hi is the 
surface field^. 

2) We have neglected the fluctuations varying in 
the lateral plane with wavelengths longer than the 
three-dimensional ^. Thus the capillary condensation 
transition has been treated at the mean-field level, 
leading to Fig. 13. In idealized conditions, there should 
be composition-dependent crossovers from the Ising 
behavior in three dimensions to that in two dimensions. 

3) Nucleation and spinodal decomposition should take 
place between plates and in porous media if T is changed 
across the capillary condensation line outside the solvent 
coexistence curvei>S.. 

4) For neutral colloids, the attractive interaction arises 
from overlap of composition deviations near the colloid 
surfaces. It is intensified if the component favored by the 
surfaces is poor in the reservoir. A bridging transition 
further takes place at lower r between strongly adsorbed 
or wetting layers of colloid particles^^i^. 

5) Our local functional theory has been used when ?/; 
varies over a wide range in strong adsorption. It can 
be used in various situations. For example, dynamics 
of colloid particles in near-critical fluids can be studied 
including the hydrodynamic flow. So far, phase sep- 
aration in near-critical fluids has been investigated in 
the scheme of the ip'^ theory with constant renormalized 
coefficients^. However, the distance to the bulk critical- 
ity (the parameter w in Sec. HI) can be inhomogeneous 
around preferential walls or around heated or cooled 
walls. 

6) We should further investigate the ion effects in 
confined multi-component fiuids. In such situations, the 
surface ionization can depend on the ambient ion den- 
sities and composition^- . A prewetting transition then 
appears even away from the solvent coexistence curve, 
where the degree of ionization is also discontinuous. We 
are interested in how the ionization fluctuations affect 
the ion-induced capillary condensation transitior^i^. 
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Appendix A: Mechanical equilibrium 

In one-dimensional situations, the zz component of the 
stress tensor due to the composition deviation is^ 

n.. - - / + C(|^f /2 - - CVlVf /2, (Al) 

where C" = dC /d4>. The mechanical equilibrium condi- 
tion is dWzz/dz — 0, so 1X22 is a constant independent of 
z. We further use Eq.(2.11) to eliminate the term pro- 
portional to ip" to obtain 



n^. - ^A^co - / + ci^Af /2. 



(A2) 



Thus — ~ /m at z = D/2. The osmotic 

pressure is given by 

n = Hzz - (V'ooMoo - /oo), (A3) 

so we find n = -dVL/dD in Eq.(2.39). 

Appendix B: Relationship to the Schofield, 
Lister, and Ho linear parametric model 

The linear parametric model**^ provides the equation of 
state and thermodynamic quantities of Ising systems in 
compact forms^i^ for detailed discussions on this model. 
It uses two parametric variables, f and 9, with f > and 
\0\ < 1; f represents a distance to the critical point and 
9 an angle around it in the ip-r plane. Here r should 
not be confused with r in Eq.(3.5). In this model, homo- 
geneous equilibrium states are supposed. The reduced 
temperature r, the magnetic field h, and the average or- 
der parameter -0 are expressed in terms of f and 9 as 



r 
h 



f(l 
Oof 
Cor' 



(Bl) 
(B2) 
(B3) 



Here ao and cq are positive nonuniversal constants, while 
he is a universal number. The case = corresponds to 
T > and /i = 0, 6* = ±l/6c to r = 0, and 6* = ±1 to the 
coexistence curve (/i = and r < 0). We may calculate 
various thermodynamic quantities from these relations in 
agreement with the asymptotic critical behavior. Though 
he is arbitrary within the model, h\ was set equal to 



This choice yields simple expressions for the critical am- 
plitude ratios in close agreement with experiments. The 
linear parametric model in Eqs.(Bl)-(B4) is exact up to 
order e^— . The two-scale-factor universality^"'''"^ further- 
more indicates that the combination (aoCo)'?o ^^e coef- 
ficients in Eqs.(B2) and (B3) should be a universal num- 
ber, where ^0 is the microscopic length in the correlation 
length ^ = \qt~" for r > and /i = 0. 

Our model in Eqs.(3.4)-(3.9) closely resemble the liner 
parametric model as regards the thermodynamics of ho- 
mogeneous states. Note that h in Eq.(B2) corresponds 
to [ilksT in Eq.(3.20) and f to w in Eq.(3.9). For our 
model, we may introduce the angle variable 9 by 



= sign(^)&oiVr^, 



(B5) 



where S is defined in Eq.(3.17). We set 9 — ±1 on the 
coexistence curve so that 



hi 



1 + 1/cr?^ 1.58, 



(B6) 



where a is given by Eq.(3.22). Then t, and 1/' are 
expressed in terms of w and 9 as 



T = w{\-hl9^), 



(B7) 



1 - Ai6i2 



(B9) 



^ = {ho/^/C2)w"9, 



where the coefficients Ai, A2^ and C3 are expressed as 

0.23, 



A2 



5a6^/[5a + (2-a)cr] 
(1 - 2P)hl = 0.55, 



(BIO) 
(Bll) 



C3 = 6^[5a+(2-a)a]Ci/[6VC2eo']- (B12) 

As differences between our model and the linear para- 
metric model, 6q in Eq.(B6) is larger than h^ in Eq.(B3) 
by 16% and there appears the extra factor. 



Z{9)^{l^Ai9^)/{l-A29^), 



(B13) 



(^-3)/[(^-l)(l- 2/3)] = 1.36. 



(B4) 



in the right hand side of Eq.(B8). Our model consid- 
erably deviates from the parametric model close to the 
coexistence curve mainly because of Z(l) = 1.73. The 
correlation length and the susceptibility on the co- 
existence curve are thus underestimated in our theory as 
in Eqs.(3.26) and (3.27). Also in our model the combi- 
nation C3(6o/"\/C2^)Co for coefficients in Eqs.(B8) and 
(B9) is a universal number from Eq.(3.11) in accord with 
the two-scale-factor universality. 
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